TEST OF THE INVERSE SQUARE LAW THROUGH 
PRECESSION OF ORBITS 



N.I. Kolosnitsyn and V.N.Melnikov 

Center for Gravitation and Fundamental Metrology, VNIIMS 
and Institute of Gravitation and Gosmology, 
Peoples' Friendship University of Russia, 
3-1, M.Ulyanovoi Sir., 119313, Moscow, Russia 
e-mail: melnikov@rgs.phys.msu.su 

rn"; 

O I Abstract 
O . 
(N ■ 

I Using precession of orbits due to non- Newtonian interaction between two celestial bod- 
■ ies and modern tracking data of satellites, planets and a pulsar we obtain new more precise 
limits on possible Yukawa-type deviations from the Newton law in planets (satellites) radii 
ranges. 

1. Non-Newtonian interactions, or range variations of G. 
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. Many modified theories of gravity and unified theories [1, 2, 3] predict some deviations 
^ from the Newton law (inverse square law, ISL) or composition- dependent violations of the 
O ' Equivalence Principle (EP) due to appearance of additional fields or new possible massive 
2 particles (partners) [5, 6]. Experimental data limit the existence of these effects at nearly 
all ranges at some level, but the most poor data are available at less than millimeter and 
also at meters and hundreds of meters ranges. The only positive result of existence of such 
, deviation in the range of 20 to 500 m was obtained by Achilli et al. using an energy storage 
plant experiment with gravimeters. They found the deviation from the Newton law with 
[ the Yukawa potential strength a between 0.13 and 0.25. This result contradicts other 
^ ' experimental data and needs to be verified in similar or other independent experiments, 
probably in space ones [6, 4, 17]. 

In the Einstein theory G is a true constant. But, if we think that G may vary with 
time [3, 5], then, from a relativistic point of view, it may vary with distance as well. In GR 
massless gravitons are mediators of the gravitational interaction, they obey second-order 
differential equations and interact with matter with a constant strength G. If any of these 
requirements is violated or we consider quantum corrections to the classical theory, we 
come in general to deviations from the Newton law with range (or to generalization of 
GR). 

In [6] several classes of such theories were analysed: 

1. Theories with massive gravitons like bimetric ones or theories with a A-term. 

2. Theories with an effective gravitational constant like the general scalar-tensor ones 

[5]. 

3. Theories with torsion. 

4. Theories with higher derivatives (4th-order equations etc.), where massive modes 
in a propagator appear leading to short-range additional forces. 

5. More elaborated theories with other mediators besides gravitons (partners), like 
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supergravity, superstrings, M-theory etc. 

6. Theories with nonhnearities induced by any known physical interactions (Born- 
Infeld etc.), where an effective mass appears. 

7. Phenomenological models, where the detailed mechanism of deviation is not known 

(fifth or other force). 

In all these theories some effective or real masses appear leading to Yukawa-type 
deviation from the Newton law, characterized by strength a and range A or to a power 
law deviation (see [6]). 

There exist some modcl-dcpcndant estimations of these forces. The most well-known 
one belongs to Scherk (1979) from supergravity where the graviton is accompanied by 
a spin-1 partner (graviphoton) leading to an additional repulsion. Other models were 
suggested by Moody and Wilczek (1984) - introduction of a pseudo-scalar particle - 
leading to an additional attraction between macro-bodies with the range 2 • 10~^ cm 
< A < 20 cm and strength a from 1 to 10~^° in this range. Another supersymmetric 
model was elaborated by Fayet (1986, 1990), where a spin-1 partner of a massive graviton 
gives an additional repulsion in the range of the order 10^ km and a of the order 10^^^. 

A scalar field to adjust A was introduced also by S. Weinberg in 1989, with a mass 
smaller than 10~^ eV/c^, or a range greater than 0.1 mm. One more variant was suggested 
by Peccei, Sola and Wetterich (1987) leading to additional attraction with a range smaller 
than 10 km. Some p-brane models also predict non-Newtonian additional interactions, in 
particular in the mm range, what is intensively discussed nowadays in the hierarchy and 
braneworld models. About PPN parameters for multidimensional models with p-branes 
see [8]. 

In this paper we consider limits on a possible Yukawa type additional interaction in 
the ranges of the order of planets (satellites) distances from the Sun (the Earth) using 
the precession method and satellites plus planets tracking data. 



2. Basic formulas. 



As it is known, in two bodies problem an orbit of a celestial body is closed for only 
two interaction potentials [9]. They are: 

(1) the Newtonian potential: C/ ~ 1/r and 

(2) U 

In other cases the orbit is not closed and a pericentcr precession is observed. In 
particular, the deviation from the Newton law due to the Yukawa additional interaction 



C/' = ^^aexp(-r/A) (1) 
r 

entails a precession of an orbit. 

In a general case the precession magnitude due to a small perturbation, described by 
a potential 5U , is equal to (see 
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Integration is done over a "non-perturbed" trajectory. Here mi is a mass of one 
celestial body, 1112 is the mass of another celestial body, M = mr'^ip is the integral of 
motion (the angular momentum), m = 17111712/ {rrii + 1712) is the reduced mass, 

5U — a{Gmim2/r) exp{—r / X). 

The "non-perturbed" trajectory is described by the expressions: 

^ = ^1- > ^=1 + ^7^ P=~F^ = a(l-e . 3 

l-|-ecos(/7 m[Gmim2) mGmim2 

where e is an eccentricity, a is a semi-major axis. 

After differentiating the right-hand side of Eq.(2) with respect to M we obtain 



Taking into account Eq. (3) we get 



dr r 



dM M 



e + cos (f 
e(l -|- ecosifi) 



(5) 



After substitution of (5) in (4) we obtain 



" ^^iP ' '''' 



e + cos (/? 



e + cos (/? 



e(l -|- e cos (^) Ay e(l -|- e cos (p) 



dip (6) 



As a result, we have 



2 / exp(— r/A) f r 
Of/? = a- 



e i (1-1- ecos(/3)2 I A 



2e -|- (1 -|- e^) cos(/3 



e -I- cos 



(7) 



where 



r a 1 — e 



A A (1 -I- ecos</?) 



(8) 



3. Estimation of effects. 
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Using Eq. (7) and proposed data on 5(p error for the LAGEOS Satellites, the inner 
planets, Moon and the binary pulser B1913+16, we calculated curves a{\), which deter- 
mine a bound in the plane (a, A) between two domains, where the Yukawa interaction (a 
new nonnewtonian force) is forbidden by experiment and where it is not. The sensitivities 
to Yukawa interactions - the united curve a{\) described by Eq.(7) - are shown in Fig.l 
as the domain 1998 for the parameter A in the range from 1 • 10^ m to 1 • 10^^ m. 

We took data for LAGEOS and LAGEOS II from [11] related to the relativistic Lense- 
Thirring precession. The data contain information on classical (Newtonian) and relativis- 
tic precession. The last one equals to ~ 57 milliarcscconds per year for LAGEOS II. The 
error for test of the Lcnsc-Thirring effect is equal to 20%. Therefore, we took the error 

= 11.4 milliarcscconds per year. The LAGEOS II eccentricity e = 0.014. Using the 
data we obtained ajnin — 1-38 • 10"^^ for A = 6.081 • 10^ m. Our result is in a very good 
agreement with the estimation for a ~ 1 • 10~^^ obtained by lorio (see [20]. 

For the Moon there are high precision data of Laser Ranging (see, for example [12]). 
In Ref. [12] geodetic precession (first calculated by de Sitter), equaled to 19 ms per year, 
was determined with the error 0.9%. The absolute error was equal to 1, 71 • 10^"^" /y. We 
used this estimation and obtained amin — 3.64 • 10~^^ for A = 1, 92 • 10^ m. 

Very precise information on the Mercury precession is contained in the article of Pitjeva 
[13]. Using radar observations of Mercury for the time span 1964-1989 and taking its 
topography Pitjeva has estimated systematic errors for the perihelion motion and a new 
random error for this motion as 0.052" /q/. From this it follows that amin = 3.57 • 10~^° 
for (A = 2.89 -10^° m). 

The most precise data on the Mars precession are obtained from observations of the 
landing spacecraft Pathfinder (see [14] ). Pitjeva has analysed joint observation data of 
Viking and Pathfinder missions and has obtained more precise estimation for precession 
of the Mars [19]. However a sensitivity to the Yukawa force in this case is very small and 
the sensitivity curve for Mercury dominates. 

We also used data for the binary pulsar PSR 1913-1-16. The pulsar has an orbit 
with the high eccentricity e = 0.617 and a big precession 4, 23° /y ± 0, 000007°/y. Using 
parameters of the pulsar in [15] we calculated amin = 6.409 • 10~^^ for A = 7.515 • 10^ m. 

Fig.l shows our data (LAGEOS, MOON, PSR1913, MERCURY) as the new exper- 
imental constraint on the coupling parameter a as a function of the parameter A taken 
from the reference [16]. 

4. Conclusions 

Here, we presented some new limits on possible deviations from the inverse square 
law using planets and satellite tracking data. They improve the existing ones by several 
orders at satellite (~ 6 • 10^ m) and planets distance from the Sun (~ 1 • 10^^ m) ranges. 
Some proposed space projects (such as SEE) make possible to improve these estimations 
to five orders at a meter range (see [17] and [10] ). Limits at less than mm ranges see in 
[18]. 
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